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Abstract. This article is concerned primarily with the 
intrinsic anisotropy in the angular distribution of 2281 
gamma-ray bursts (GRBs) collected in Current BATSE 
Gamma-Ray Burst Catalog until the end of year 1998, 
and, second, with intrinsic anisotropies of three subclasses 
("short", "intermediate", "long") of GRBs. Testing based 
on spherical harmonics of each class, in equatorial coor- 
dinates, is presented. Because the sky exposure function 
of BATSE instrument is not dependent on the right as- 
cension a, any non-zero spherical harmonic proportional 
either to P™ (sin 5) sin ma or to P,™ (sin S) cos ma with 
m 7^ (d is the declination), immediately indicates an 
intrinsic non-zero term. It is a somewhat surprising re- 
sult that the "intermediate" subclass shows an intrinsic 
anisotropy at the 97% significance level caused by the high 
non-zero P 3 X (sin S) sin a harmonic. The remaining two sub- 
classes, and the full sample of GRBs, remain isotropic. 
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- gamma rays: bursts 



1. Introduction 

The discovery of an anisotropy in the angular distribution 
of 2025 gamma-ray bursts (GRBs) collected in Current 
BATSE Gamma-Ray Burst Catalog ( |Meegan et al. 1998| ) 
up to the end of year 1997, has recently been an- 
nounced ( Balazs ct al. 199"§| ). Of course, the existence 
of an observed anisotropy is expected from the non- 
uniform sky exposure function of BATSE instrument 
(|Fishman ct al. 1994|). Nevertheless, Balazs et al. (1998) 



argue that the detected anisotropy should not be a pure 
instrumental effect; i.e. there also exists an intrinsic 
anisotropy (see also Balazs et al. 1999). The exact details 
of the nature of the intrinsic anisotropy remains open to 
debate. 



Recently, Horvath (1998) and Mukherjee et al. (1998) 
have proposed that the long GRBs should be further sep- 
arated into two classes; the limit separating the "long" 
subclass should be at Tgo = 10 s, where Tgo is the "90%" 
duration (Fishman et al. 1994). Hence, in the studies of 
subclasses it is required to consider three subclasses sepa- 
rately ("short" subclass with Tgo < 2 s; "intermediate" 
subclass with 2 s < Tgo < 10 s; "long" subclass with 
T 90 > 10 s). 

In principle, for GRBs, there are 16 possible types 
which can occur, because for " all" GRBs, and for the three 
subclasses as well, there are two eventualities: either they 
are distributed isotropically or they are not. 

The purpose of this paper is to clarify the situation 
concerning these possibilities. We will again use the anal- 
ysis based on the spherical harmonics (Balazs et al. 1998). 

In the study by Balazs et al. (1998), an analysis of 
GRBs having a peak flux greater than 0.65 photons/ (cm 2 
s) was performed. Hence, the GRBs having a peak flux 
smaller than 0.65 photons/(cm 2 s) were omitted from the 
sample. This truncation did not give any new result; ex- 
cept for the fact that the number of GRBs was lowered, 
and hence the corresponding significance levels were also 
lowered. The High-Energy versus Non-High-Energy sepa- 



ration (Pendleton et al. 1997) did not give any new result 
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either; hence, we will not use this separation here. 

It is well-known that there is no general agreement con- 
cerning the significance level required to reject a null hy- 
pothesis. In an ad hoc manner, we will take, in accordance 
with the propositions of standard statistical textbooks 
( [Trumplcr fc Weaver 1953| , [Kendall fc Stuart 1969[ ), for 
the "doubtless" or "definite" rejection, the > 99.9% sig- 
nificance level. (This level of significance is obtained for 
all GRBs in the sample examined by Balazs et al. 1998; 
but together with the instrumental effects.) We will treat 
a significance level of > 99% as a practically sure result, 
with the > 95% significance level result being considered 
as noteworthy. 
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2. Mathematical considerations 

The key ideas of the analysis by Balazs et al. (1998) were 
based on spherical harmonics up to the quadrupole terms. 
In this section, we generalize these considerations to in- 
clude higher order harmonics. In addition, we introduce 
new statistical tests. 

As we note in the introduction, it is necessary to elim- 
inate the BATSE instrumental effects coming from the 
non-uniform sky exposure. Several methods exist. For ex- 
ample, in Briggs et al. (1996), Monte Carlo simulations 
are used: N points are randomly scattered on the sky 
many times, and the scattering is "deformed" in accor- 
dance with the sky exposure function; N is the number 
of GRBs. Then the spherical harmonics calculated from 
these simulations are compared with the spherical har- 
monics of the observed BATSE data. Moreover, in Briggs 
et al. (1996) a further quite simple method is proposed: 
to use the equatorial coordinates a and 5 (0 < a < 2tt; 
—tt/2 < 5 < tt/2) instead of the Galactic coordinates, be- 
cause the sky exposure function is not dependent on a. In 
this paper we will use these coordinates and the proposed 
"equatorial" method. 

The key ideas of this method may be seen as follows. 
Let u>(a, 5) cos 5 dS da be the probability of finding a 
GRB in the infinitesimal solid angle cos 6 dS da. If we 
assume that there is an isotropy in the sky distribution 
of GRBs, then ui sh ould reflect the sky exposure function 
( Balazs et al. 199"s| ). The function u> can be decomposed 
into its spherical harmonics. One then has 



EE^ 

n — 1 m— 1 



u>(a,5) = y~V„, P n (sin 5)- 



(sin5)(w n! _ m sin ma + w„ jra cos ma), 



(1) 



where P™(sino") are the Legendre polynomials. 

If there is an intrinsic isotropy, then uj(a,5) = w(5), 
and one must have w ni _ m = uj n .m = for any n > m > 1. 
Conversely, if even one single w ni ± m for m > 1 were non- 
zero, then there would be an intrinsic anisotropy. This is 
the key idea of this method. Simply, if to is not dependent 
on a, then all m =/= harmonics must be identically zero 
independent of the concrete form of function uj(5). Then, 
of course, this theoretical prediction should be compared 
to the observational data. 

Concerning the m = terms one ma y use either the 
Monte Carlo method ( Briggs et al. 1996 ), or - if possible 
- to calculate from the known w(<5), numerically or even 
analytically, the expected theoretical values of co n ,o- They 
will not be zeros in the general case. (They would be zeros, 
if lu were constant; i.e. if the sky exposure function were 
constant.) 

In this paper we will not study these terms in detail, 
and an eventual observed non-zero Lu n fl will simply be 



assumed to be that caused exclusively by the instrumental 
effects. 

In our case the null hypothesis of intrinsic isotropy will 
be identical to the assumption a;(a,<5) = lo(S); i.e. to the 
assumption that for all w n ,± m = for m ^ 0. Hence, to 
test the null hypothesis of intrinsic isotropy one has to test 
these theoretically expected zeros for the given spherical 
harmonics with m/0. 

This may be done similarly to Balazs et al. (1998). This 
means that, if one has N GRBs with measured [ctj, 5j] 
(j = 1,2, N) coordinates, one has to calculate for the 
given n,m the values (Peebles 1980, Eq. (46.1); Balazs et 
al. 1998) 



N 



^2 p n (sin Sj) cos ma j 
and 



\ 



(P™ (sin Sj ) cos may ) 5 



(2) 



(3) 



i=i 



(For n, —m one simply has to substitute cosine by sine.) 
Then one can apply Student's t test, where the expected 
mean is zero. Here t is simply the ratio of quantities from 
Eqs. (2-3), because for N » 1 (this is always the case 
here) one may use N — 1 ~ N. Note also that there are 
no problems which will arise from the normalization con- 
stant of a spherical harmonic. (If, instead of P™, one uses 
const. P™, where const, is an arbitrary positive number, 
then t will not change, if N — 1 ~ JV, and if the expected 
mean is zero.) Technically, there can also occur several 
computational numerical problems with the definition of 
P™ itself for large n on the computer (see Press et al. 
1992, Chapt. 6.8). In our calculations, up to the used n, 
we encountered no such problems. We used the algorithm 
described in Chapt. 6.8 of Press et al. (1992). 

We will also use a modification of this test. 
From the symmetry properties of P™(sin S) cos ma and 
P™ (sin 8) sin ma it follows that for m ^ in exactly one 
half of the sky the sign of the harmonic is positive, and 
in the other half it is negative. For m = this is true for 
odd n only, but not for even n. Hence, except for even 
n with m = 0, one may proceed as follows. One simply 
calculates the number fc Q b s of GRBs having - say - neg- 
ative values. Then one may apply the Bernoulli test de- 
scribed in Meszaros (1997) and Balazs et al. (1998). (It is 
essential to note that the half of the sky, where the given 
harmonic has the given sign, need not be a "connected" 



compact region (Balazs et al. 1998).) Then for the given 
kobsobs the calculation of the significance level is described 
in Balazs et al. (1998). In essence, one can take - as the 
theoretical expectation - a normal distribution with mean 
A72 and dispersion y/N/2. Hence, if |fc obs - N/2\ > \/N, 
t hen the significan ce level is larger than 2a, i.e. 95.4% 
(Balazs et al. 1998). Note that - trivially - in this "sign 



A. Meszaros et al.: On the existence of the intrinsic anisotropies in the angular distributions 



3 



test" the choice sign does not matter. Obviously, the "sign- 
test" and Student's t test need not give the same result 
for a given harmonic. 

This method based on spherical harmonics seems to be 
extremely powerful because, at least in principle, one sin- 
gle non-zero harmonic for ra/0 should be enough to re- 
ject the null hypothesis of intrinsic isotropy (Peebles 1980, 
Chapt. 46; Tegmark et al. 1996). This follows from the or- 
thogonality of spherical functions (Peebles 1980, Chapt. 
46). On the other hand, it is not excluded that - testing 
several harmonics simultaneously with two different tests 
-, some of these tests can give a higher than a% signif- 
icance level "by chance". For example, if a = 95, then 
one may expect that 5% of tests can give a greater than 
95% significance level "by chance". Therefore, to be max- 
imally careful, we will reject the null hypothesis only in 
the case, when the number of tests giving > a% signifi- 
cance level is far above (100 — a)% fraction. Concretely, 
this means the following. We take the absolute value of 
the difference between the number of tests giving > a% 
significance level and the theoretically expected number 
of tests giving > a% significance level. This value must 
be so large that the probability of a chance of this differ- 
ence is smaller than (100 — a)%. This probability is simply 
calculated from the binomial (Bernoulli) distribution. 

In our case we have infinite harmonics and hence in 
principle we may apply an infinite number of tests. Never- 
theless, due to the discrete character of measured oj(a, 5), 
there is a limitation to the harmonics, which can be tested. 
Having the angular coordinates of N GRBs, we have 2N 
independent measured quantities. Hence, clearly, we can- 
not obtain more than 2N independent harmonics. Which 
harmonics are these? In addition, even for these 2N inde- 
pendent harmonics there may occur a further problem. Of 
course, the mean, variation and the t for a given harmonic 
can be calculated from the positions using Eqs. (2-3). But 
what is the accuracy of this harmonic calculated for finite 
N compared with the ideal case, when one would have 
N — > oo? In other words, is the calculated t a "true" es- 
timator of the Student t? Obviously, we would need to 
compare the measured t obtained for the "infinite" case 
with the theoretically predicted zero mean. 

Several authors have investigated the answers to 



these and similar questions ( 


Irumpler & Weaver 1953, 


Kendall & Stuart 196S|, 


Bracewell 1978 




Peebles 198C, 


Press et al. 1992 




Tegmark et al. 1996 


). Note that some 



conclusions can be deduced for this case immediately 
from the well-discussed case of discrete Fourier transform 
(cf. Bracewell 1978, Chapt. 14; Press et al. 1992, Chapts. 
12-13); in a we have in fact a Fourier decomposition. 
The key result of all these studies is the statement 
that for the lowest-order harmonics (dipole, n = 1; 
quadrupole, n — 2; etc..) the Student t calculated 
from the obser vational data can be considered as if i t 
oo jHorack et al. 1993| , flegmark et al. 1996|) . 



were N — > 
Hence, for 



these low-order harmonics Student's t test 



gives the correct conclusion. This result also suggests that 
the 2N independent harmonics, which should be tested, 
should be the lowest possible ones with n < (y/2N — 1). 
(Clearly, up to n, we have 1 + 3 + ... + (2n + 1) = (n + l) 2 
harmonics.) In addition, one should note that, even for 
n > \/2N — 1 some harmonics, being expressed by the 
lower ones, can also be used with care. For example, in 
Tegmark et al. (1996), N was = 1122 and even n = 65 was 
discussed. Nevertheless, in this paper, in order to avoid 
any complications with higher n, we restrict ourselves to 
n < (V2N- 1). 

A rough estimate of the inaccuracy of quantities in 
Eqs. (2-3) for the given n,m,N may be quite simply 
performed. The typical angular scale belonging to n is 
~ ir/n (Peebles 1980, Chapt. 46), and hence the corre- 
sponding solid angle is ~ {n/n) 2 steradians. Clearly, on a 
solid angle ~ Ait/N steradian the discreteness is crucial, 
and hence for {tt/u) 2 ~ (Att/N) the inaccuracies can be 
even of order unity (they need not be, but they can be 
for the worst case). All this suggests that for n the in- 
accuracy coming from discreteness should be maximally 
~ (47T / TV) / (tt I 'n) 2 ~ n 2 /N. Of course, in the inaccuracy 
formula ~ n 2 /N there is an uncertainty concerning the nu- 
merical factor itself. (We have roughly taken 4 ~ 7r; in Eq. 
(3) cos 2 na may be expressed by cos 2na, and hence here 
one will have generally a four-times larger uncertainty; 
the typical ir/n angle belonging to n is a rough estimate; 
etc.) Recently, this question has usually been studied by 
numerical simulations (cf. Tegmark et al. 1996, and ref- 
erences therein). But, classical mathematical treatments 
are also known (cf. Boas 1983, Chapt. 16.6; Press et al. 
1992, Chapt. 13.4; Bracewell 1978, Chapt. 14). All these 
studies suggest that for n <C \/2N the discreteness leads to 
negligible inaccuracies, which is in accordance with our es- 
timate ~ n 2 /N. On the other hand, Tegmark et al. (1996) 
imply that even for n 2 ~ N the analysis based on spherical 
harmonics can be good. This suggests that the inaccuracy 
should be much smaller that n 2 /N. (Note that, trivially, 
the dependence on m need not be mentioned especially. 
Clearly m < n, and there are 2n + l harmonics for a given 
n. Then the inaccuracy for a given n will simply be the 
largest inaccuracy among these 2n + 1 harmonics.) We 
will not go into the details of this discussion, and in this 
paper we will simply assume that the inaccuracy is maxi- 
mally ~ n 2 /N. We will also keep in mind the fact that the 
estimate of this inaccuracy is never complete. Therefore, 
one should be careful with any conclusions derived, for in- 
stance, when n 2 -C TV does not apply. To avoid problems 
with higher harmonics, we will carry out a further dras- 
tic truncation. We will, ad hoc, consider only n < y/N /3. 
This means that instead of the possible 2N harmonics, we 
will study only ~ N/9. This (together with the avoiding 
of m = terms) means that we will test only ~ (5 — 6)% 
of all allowed independent harmonics. 

Concerning the "sign-test" the accuracy seems to be 
much better than for Student's t test. By this we mean 
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to say that there is no problem with the correctness of a 
conclusion coming from the Bernoulli test. Therefore, the 
" sign-test" may be used for the same harmonics for which 
Student's t test was used. 

For the used values of n the positional errors (of sizes 
~ 1 — 4 degrees, Fishman et al. 1994) should give no fur- 
ther complications. We will consider only such values of n, 
when the typical angular size belonging to an n (~ ir/n; 
Peebles 1980, Chapt. 46) will be much larger than the size 
of the positional error. 

The detailed discussion presented by Tegmark et al. 
(1996) shows that for such values of n no complications 
should arise from the positional errors. 

3. Results 

In order to test the isotropy of 2281 GRBs col- 
l ected in Current B ATSE Gamma-Ray Burst Catalog 
( Mccgan et al. 199^ ) until the end of 1998, we calculated 
their spherical harmonics up to n = 15. This maximum 
n value was chosen in accordance with the discussion in 
Sect. 2 , where we restricted ourselves to n < y/~N /3. Here 
N = 2281, and s/N/3 = 15.9. In other words, instead 
of 2 x 2281 = 4562 possible independent spherical har- 
monics we tested only (16 2 — 16) = 240 both with Stu- 
dent's t test and sign tests. The results are collected in 
Table 1, in which we present all dipole and quadrupole 
terms (in order to compare them with Balazs et al. 1998), 
but for larger harmonics only the terms, for which either 
|t| > 1.96 or fc obs < 1093 or k oha > 1188. For |f| > 1.96 
there is a smaller than 5% probability that there is still an 
isotropy; for \t\ > 2.58 this probability is smaller than 1% 
(Trumplcr & Weaver 1953). For the "sign test" one needs 
\k bsobs — N/2\ > y/ N to have a > 95.4% significance level 
for anisotropy, where fcobs is the observed number of GRBs 
having the given harmonic negative value. We also present 
the to = cases for comparison as well, but these values 
will not be taken as indications of intrinsic anisotropy. For 
m = the "sign test" is applied only for odd n, of course. 

Table 1 shows that there are three to = components 
with |*| > 1.96 (n = 1,2,12). The biggest value is the 
n = 2, to = quadrupole term with t — 4.26. This means 
that there exists an anisotropy with a certainty; the proba- 
bility that this quadrupole term is non-zero due to chance 
is much smaller than 0.1% probability (* = 3.29 corre- 
s ponds to a 0.1% p robability). In fact, this is nothing new 
( Balazs et al. 199*^ ); it again confirms that the observed 
distribution of all GRBs is without doubt anisotropic. 

We will interpret these three terms as the anisotropy 
that follows exclusively from the non-uniform sky exposure 
function of BATSE instrument. 

There are no further dipole and quadrupole 
anisotropies above |*| > 1.96; the sign test also gives 
no such anisotropies. On the other hand, there are 19 
m 7^ harmonics above the 95% significance level. (In 12 
cases |*| > 1.96; in 7 cases the sign test gives anisotropy 



Table 1. Results of Student's * and the sign tests, re- 
spectively, of 2281 GRBs from spherical harmonics up to 
n = 15. All dipole + quadrupole terms are presented, but 
for higher order harmonics only the components with ei- 
ther |*| > 1.96 or fc b s < 1093 or fc obs > 1188. * (**) 
means that the test gives a > 95% (> 99%) probability of 
anisotropy. For to > (to < 0) one has the cosine (sine) 
term in Eq. (1). For n even with to = the sign test is 
not applied. 



n 


rn 


* 


fcobs 


1 





*2.29 


1098 


1 


1 


—0.44 


1155 


1 


— 1 


0.12 


1148 


2 





* * 4.26 




2 


1 


—0.88 


1149 


2 


— 1 


—0.30 


1149 


z 


o 
1 


— u.zy 


1 1 

llo i 


2 


-2 


0.85 


1137 


4 


4 


0.89 


*1092 


5 


4 


*2.20 


1094 


9 


-6 


0.53 


*1092 


9 


8 


* * 2.67 


* * 1072 


9 


-8 


*2.44 


1106 


10 


-1 


1.42 


*1089 


10 


-10 


* - 2.14 


1184 


11 


8 


* * 2.88 


1136 


11 


-8 


*2.03 


1095 


12 





*2.31 




12 


-12 


* * 2.69 


1094 


13 


8 


* * 2.81 


* * 1072 


13 


9 


*2.18 


1147 


14 


-9 


0.09 


*1092 


14 


-12 


* * 3.36 


1103 


15 


8 


*2.39 


1115 


15 


9 


*2.43 


1110 


15 


11 


0.28 


*1091 



above the 95% significance level; for n = 9, to = 8 and 
n = 13, to = 8, respectively, both tests give anisotropy.) 
In addition, for five harmonics one has |*| > 2.58. For two 
of them one has fc bs < 1093 as well. Concerning the 95% 
significance level, one expects that in 2 x 240/20 = 24 
cases, one can obtain this "by chance", too. This is more 
than the 19 cases obtained, and the difference could well 
be due to chance. Concerning the 99% significance level, 
one expects 4.8 cases theoretically, and we obtain 7 ones. 
This can also be due to chance. The dispersion arising 
from the binomial distribution is \/4.8 x 0.99 = 2.18; this 
is practically identical to (7 — 4.8) = 2.2. There is also 
a * = 3.36 value, which should give a higher than 99.9% 
significance level. This single value is hardly enough to 
reject the null hypothesis. One expects 0.48 such cases. 
The dispersion is 0.69; hence one single test may well 
occur "by chance" as well. Hence, we conclude that the 
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Table 2. Results of the Student t + sign tests, respec- 
tively, for 419 short GRBs from spherical harmonics up to 
n = 6. The notation is the same as in Table 1; for k Q k s 
here we needed either fc Q b s < 190 or fc Q b s > 231. 



n 


m 


t 


^obs 


1 





0.38 


207 


1 


1 


-0.61 


211 


1 


-1 


1.52 


200 


2 





* * 2.69 




2 


1 


-1.07 


218 


2 


-1 


-0.05 


207 


2 


2 


-0.08 


201 


2 


-2 


1.74 


199 


3 


1 


0.98 


*186 


5 


-1 


1.60 


*188 


5 


2 


0.78 


*189 


5 


3 


* - 2.19 


229 



null-hypothesis of intrinsic isotropy for all GRBs should 
not be rejected. 

In Current BATSE Catalog there are 1670 GRBs hav- 
ing measurements both for the positions and T90. (The du- 
rations are catalogued for GRBs detected before August 
1996.) Among these, 419 objects have T 90 < 2 s ("short" 
subclass), 253 objects have 2 s < T 90 < 10 s ("interme- 
diate" subclass), and 998 objects T 90 > 10 s ("long" sub- 
class). For these three subclasses we repeated the above 
procedure. 

The results for the 419 short GRBs having T 90 < 2 s 
are collected in Table 2. (Here N = 419, and therefore we 
consider only n < 6.) Again there is a clear non-zero term 
for n = 2, m = expected from the BATSE 's sky exposure 
function. There are no dipole and quadrupole terms above 
\t\ > 1.96. On the other hand, there arc 4 further tests 
determining non-zero harmonics above the > 95% signifi- 
cance level. Nevertheless, from the 84 tests performed here 
one expects theoretically that 4.2 tests will be above the 
95% significance level "by chance". Again, similarly to all 
GRBs, the null-hypothesis of intrinsic isotropy of short 
GRBs should not be rejected. 

The results for the 253 intermediate GRBs having 2 s 
< T 90 < 10 s are collected in Table 3. (Here N = 253, and 
therefore we consider only n < 5.) Interestingly, there are 
no non-zero m = terms expected from the BATSE sky 
exposure function. There are 6 tests determining non-zero 
harmonics above the > 95% significance level. From the 60 
tests done here one expects theoretically that 3 tests will 
be above the 95% significance level "by chance". The ob- 
tained 6 cases - instead of the expected 3 ones - can just 
still be by chance. The dispersion is V3 x 0.95 = 1.69. 
6 - 3 = 3 is still smaller than 2 x 1.69 = 3.38. Hence 
the obtained 6 tests instead of the expected 3 can still be 
by chance with a probability larger than 5%. Neverthe- 



Table 3. Results of the Student t + sign tests, respec- 
tively, for 253 intermediate GRBs from spherical harmon- 
ics up to n = 5. The notation is the same as in Table 1; 
for fc bs here we needed cither fc b s < HI or fc b s > 142. 



n 


m 


t 


^obs 


1 





1.33 


116 


1 


1 


-0.83 


136 


1 


-1 


* - 2.24 


142 


2 





1.74 




2 


1 


1.17 


120 


2 


-1 


-0.68 


132 


2 


2 


-0.49 


129 


2 


-2 


-0.43 


132 


3 


-1 


* * —3.37 


* * 154 


5 


-1 


* - 2.25 


130 


5 


1 


0.20 


*108 


5 


4 


*2.25 


112 



less, we find that for the intermediate subclass the null 
hypothesis of isotropy can be rejected at the > 95% sig- 
nificance level. This follows from the large value of t and 
/c bs for n = 3, m = — 1. Student's t test yields a proba- 
bility of zero value for this harmonic of 0.08%. From the 
"sign-test" this value is even smaller, because - from the 
Bernoulli distribution - the zero value is rejected at the 
(154 — 126.5)/7.95 = 3.46c significance level; the proba- 
bility of a chance is only 0.06%. From the considered 60 
tests, we obtain 2 cases giving a > 99.92% significance 
level; the expected value is maximally 0.48. Then the dis- 
persion is 0.7, and (1-0.48) = 1.52 > 2x0.7. Hence, there 
is a 1.52/0.7 = 2.17er - i.e. 97% - probability that this is 
not by chance. Taking into account this estimate we con- 
clude that this single spherical harmonic alone is enough 
to reject the assumption of isotropy at a higher than 95% 
significance level. This significance level is minimally 97%. 

The results for the 998 GRBs having T 90 > 10 s are 
collected in Table 4. (Here N = 998, and therefore here 
we consider only n < 10.) First, from these results it fol- 
lows that there are the n = 4, 7, m — anisotropy terms 
expected from the BATSE sky exposure function. There 
are no further dipole and quadrupole terms suggesting 
anisotropy on these scales. 

Second, there are 16 tests giving anisotropies between 
95% and 99% significance levels. Third, there are no tests 
giving anisotropy above the 99% probability level. Here, 
except for the m — terms, 220 test were done, and from 
them theoretically 220/20 = 11 tests should give > 95% 
probabilities by "chance". Instead of this, there are 16 
such cases, which can still be by chance. The dispersion is 
V0.95 x 11 = 3.23, and 2 x 3.23 = 6.46 > (16 - 11) = 5. 
The probability to obtain 16 cases instead of 11 by chance 
is bigger than 5%. Again, similarly to all GRBs and to the 



6 



A. Meszaros et al.: On the existence of the intrinsic anisotropics in the angular distributions 



Table 4. Results of the Student t and sign tests, respec- 
tively, for 998 long GRBs from spherical harmonics up to 
n = 10. The notation is the same as in Table 1; for k i, s 
here we needed either fc D b s < 468 or k Q b s > 530. 
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short subclass, the null-hypothesis of intrinsic isotropy of 
long GRBs should not be rejected. 

4. Conclusions 

The results of paper may be summarized as follows. 

First, the n = 2, m — spherical harmonic shows 
that there is a clear anisotropy on the significance level 
> 99.9% in the distribution of all 2281 GRBs. This fact is 
expected from the BATSE sky exposure function, and is 
interpreted as an artificial "instrumental" anisotropy. 

Second, there is a clear intrinsic anisotropy of 253 "in- 
termediate" GRBs at the > 97% significance level due to 
the n = 3 m = — 1 term. 

Third, both the 419 short GRBs and the 998 long 
GRBs, respectively, and also all 2281 GRBs do not ex- 
hibit anisotropics on statistically high enough significance 
levels. 

All these results are interesting, because the depar- 
ture from intrinsic isotropy just for the new "interme- 
diate" subclass having the smallest number of GRBs 
is surprising. Of course, the significance level should 
still be increased, because the 99% level (or even the 
99.9%) is desirable. However, the 97% significance level 
is already remarkable, and hence surely should be an- 



nounced. The refinement of significance level may fol- 
low either from the methods used in this paper and 
in Balazs et al. (1998) (for example, a better estima- 
tion of inaccuracy will not need the drastic truncation 
n < \Zn/3), or from wholly different statistical methods 
(|Kcndall fc Stuart 1969|, [Bagoly et al. 1998|). 
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